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By  performing  two  parallel  numerical  experiments  —  solving  the  dynamical  Hamiltonian  equations 
and  solving  the  Hasselmann  kinetic  equation  —  we  examined  the  applicability  of  the  theory  of  weak 
turbulence  to  the  description  of  the  time  evolution  of  an  ensemble  of  free  surface  waves  (a  swell)  on 
deep  water.  We  observed  qualitative  coincidence  of  the  results. 

To  achieve  quantitative  coincidence,  we  augmented  the  kinetic  equation  by  an  empirical  dissipation 
term  modelling  the  strongly  nonlinear  process  of  white-capping.  Fitting  the  two  experiments,  we 
determined  the  dissipation  function  due  to  wave  breaking  and  found  that  it  depends  very  sharply  on 
the  parameter  of  nonlinearity  (the  surface  steepness).  The  onset  of  white-capping  can  be  compared 
to  a  second-order  phase  transition.  This  result  corroborates  with  experimental  observations  by 
Banner,  Babanin,  Young  [13]. 


INTRODUCTION. 

Wave  turbulence  is  realized  in  plasmas,  liquid  helium, 
magnetohydrodynamics,  nonlinear  optics,  etc.  A  perfect 
example  of  wave  turbulence  is  a  wind-driven  sea.  The 
major  conceptual  difference  between  wave  turbulence 
and  “classical”  turbulence  in  an  incompressible  fluid  is 
the  presence  of  a  characteristic  dimensionless  parame¬ 
ter  pi,  characterizing  the  level  of  nonlinearity.  Turbu¬ 
lence  is  considered  to  be  “weak”  if  p  <C  1,  otherwise  it  is 
“strong”.  In  classical  hydrodynamic  turbulence  p  — >  oo. 

A  more  rigorous  definition  of  weak  turbulence  is  the 
following:  this  is  the  turbulence  which  is  well  described 
by  the  kinetic  equation  for  waves.  These  equations  are 
the  quantum  kinetic  equations  for  bosons  in  the  limit 
of  very  high  occupation  numbers.  They  were  derived  in 
statistical  physics  in  the  late  twenties  IS,  i  and  redis¬ 
covered  in  nonlinear  wave  dynamics  in  the  sixties.  The 
kinetic  equation,  describing  four-wave  resonant  interac¬ 
tion  of  gravity  waves,  was  named  after  K.  Hasselmann, 
who  derived  it  in  1962-1963  j3). 

The  theory  of  weak  turbulence  is  well-developed  Q. 
The  kinetic  equation  has  rich  families  of  Kolmogorov- 
Zakharov  (KZ)  and  self-similar  solutions,  which  can  be 
efficiently  used  for  explaining  a  wide  range  of  experimen¬ 
tal  data  Si-  However,  today  we  have  a  clear  under¬ 
standing  of  the  following  fact:  even  for  small  values  of 
p,  the  theory  of  weak  turbulence  may  be  incomplete.  In 
many  important  physical  situations  weak  and  strong  tur¬ 
bulence  coexist. 

Even  if  the  weak  turbulent  resonant  interaction  effects 
dominate  in  the  greater  part  of  space,  strongly  nonlin¬ 
ear  effects  could  appear  as  rare  localized  coherent  events. 


If  they  are  smooth  and  regular,  they  are  solitons,  qua- 
sisolitons  or  vortices.  However,  they  could  be  catas¬ 
trophic,  in  which  case  they  are  wave  collapses,  similar 
to  self-focusing  in  nonlinear  optics  or  Lagmuir  collapses 
in  plasma.  Even  rare  sporadic  collapse  events  can  essen¬ 
tially  affect  the  physical  picture  of  wave  turbulence. 

There  are  two  main  types  of  wave  collapse  events  in 
a  wind-driven  sea.  The  first  is  the  formation  of  freak 
waves;  this  is  not  a  subject  of  our  study.  The  second, 
which  is  much  more  common,  is  wave-breaking  or  white¬ 
capping,  which  is  an  essential  mechanism  of  wave  energy 
dissipation.  It  would  be  hopeless  to  develop  an  efficient 
operational  model  of  wave  forecasting  without  an  under¬ 
standing  and  a  proper  parametrization  of  this  fundamen¬ 
tal  effect.  Meanwhile,  a  reliable  analytical  theory  of  this 
phenomenon  is  still  not  developed,  while  field  and  labo¬ 
ratory  experimental  data  are  scarce.  The  most  promising 
approach  to  resolving  this  problem  is  a  massive  numerical 
experiment. 

The  most  informative  experiment  would  be  one  that 
could  provide  a  direct  numerical  solution  of  the  primi¬ 
tive  dynamic  equations  describing  the  wave  ensemble.  In 
1992,  Dyachenko,  Puslrkarev,  Newell  and  Zakharov  nu¬ 
merically  solved  2-D  focusing  NLSE  and  observed  the 
coexistence  of  self-focusing  collapses  with  weak  turbu¬ 
lence  Q.  Later  on,  the  1-D  MMT  (Maida,  McLaughlin 
and  Tabak)  model  and  its  generalizations  were  solved 
numerically  by  different  authors  (see  summary  in  [§]). 
Again,  the  coexistence  of  wave  collapses  and  weak  tur¬ 
bulence  was  verified.  In  our  article,  we  present  (as  we 
hope,  for  the  first  time  in  the  literature)  the  results  of  a 
far  more  detailed  experiment.  We  performed  the  numer¬ 
ical  simulation  of  the  evolution  of  an  ocean  swell  using 
two  different  approaches. 
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In  the  first,  we  solved  the  Euler  equations  for  the  3-D 
potential  flow  of  an  ideal  incompressible  fluid  with  a  free 
surface  in  the  presence  of  gravity.  We  used  the  Hamilto¬ 
nian  form  of  these  equations  SEl.  For  gravity  waves, 
the  parameter  of  nonlinearity  is  the  average  steepness  y. 
We  expanded  the  Hamiltonian  in  powers  of  y  up  to  order 
y4.  In  the  second  experiment,  we  solved  the  Hasselman 
kinetic  equation. 

The  comparison  of  the  results  demonstrates  qualita¬ 
tive  accordance.  Both  experiments  describe  expected  ef¬ 
fects,  such  as  the  downshift  of  the  spectrum  peak,  the 
angular  spreading  of  the  spectrum  and  the  formation  of 

Zakharov-Filonenko  spectral  tails  Fu  ~  a;-4  sum  .  To 
obtain  quantitative  coincidence  of  the  results,  we  have 
to  augment  the  Hasselmann  equation  by  an  empirical 
dissipation  term  Sdiss ,  modelling  white-capping  effects. 
We  tried  several  versions  of  this  term.  The  versions 
of  Sdiss  used  in  the  industrial  wave-predicting  models 
WAM3  and  WAM4  essentially  overestimate  the  dissipa¬ 
tion  for  a  moderate  steepness.  The  comparison  with 
dynamical  computations  shows  that  white-capping  dissi¬ 
pation  decreases  dramatically  with  decreasing  steepness 
and  that  it  is  probably  a  threshold  phenomenon,  simi¬ 
lar  to  a  second-order  phase  transition.  Similar  results 
were  earlier  obtained  in  the  field  experiment  by  Banner, 

Babanin  and  Young  E3- 


DYNAMICAL  MODEL 


In  accordance  with  recent  results  [25|,  the  dissipation 
term  should  be  included  in  both  equations. 

The  distribution  of  the  wave  action  is  described  by  the 
function  n(k,t)  =  where 


(3) 


are  complex  normal  variables.  Here  w*,  =  Vgk- 

As  the  initial  condition,  we  used  a  Gaussian-shaped 
distribution  in  the  Fourier  space: 


|a^|  =  Ai  exp 


\H\  =  10~ 


V 


k  —  kn 


<  2 A, 


k  —  k0 

Ai  =  0.92  x  10-6,  A  =  60, 
k0  =  (0;  300),  w0  =  \fgko- 


(4) 


The  initial  phases  of  all  harmonics  were  random.  The 
average  steepness  of  this  initial  condition,  defined  as  /i  = 
^2(1  Vr?|2),  was  y  ~  0.176. 

The  period  of  the  most  intensive  wave  was  To  = 
27t/V300  =  0.362.  Calculations  continued  until  t  = 
3378 T0.  We  observed  an  angular  spreading  of  the  ini¬ 
tial  spectral  distribution  together  with  a  downshift  of 
the  spectral  peak.  Level-lines  of  the  initial  and  the  fi¬ 
nal  spectra  are  presented  on  Figs.  [U  [2]  We  observed 


In  this  part  of  our  experiment,  the  surface  of  the  liquid 
is  described  by  two  functions  of  the  horizontal  variables 
x,  y  and  the  time  t:  the  surface  elevation  r](x,  y,  t )  and 
the  velocity  potential  on  the  surface  i/j(x,y,t).  In  our 
approximation,  they  satisfy  the  following  equations  Q : 


f]  =  kip  —  (V(^V^)))  —  k[r)kil>]  + 
+k(r]k[r]k'ilj ])  +  \ \r2[ri2k'tp]+ 
4fc[?72V2V>]  +  F~1['YkVk]i 


tj;  = 


-9V  -  2 


(Vi/02  —  (fc'0)" 


-[ki/)]k[T}kip]  -  [r]kip]V2il>  +  F 


(1) 


Here  k  is  the  linear  integral  operator  k  =  \AV2,  T-1 
corresponds  to  the  inverse  Fourier  transform. 

Equations  ©  are  nowadays  widely  used  in  numerical 
experiments  and  are  solved  by  different  versions  of  the 


spectral  code  [13,  DJ,  OJ,  IM  llH,  M  UJ,  l2fl  M,  \M,  123,  [24 1 . 
In  the  present  experiment,  we  solved  the  equations  in 
the  real  space  domain  27t  x  27t  using  the  finest  currently 
possible  rectangular  grid  512  x  4096,  putting  g  =  1.  The 
dissipative  terms  F~1['ykVk\  and  -F_1[7fcV’fc]  are  taken  in 
the  form  of  pseudo- viscous  high  frequency  damping.  We 
put 


r  o,  fc  <  kd, 

lk  1  ~  kd)2,k  >  kd, 

kd  =  1024,7=  5.65  x  10“3. 


FIG.  1:  Initial  spectrum  |ag|2.  t  =  0. 


Ky 

FIG.  2:  Final  spectrum  |aj|2.  i  ~  2809Tq. 


(2) 


the  following  indications  of  wave-turbulent  behavior: 
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1.  The  statistics  of  energy-capacity  spectral  modes  is 
close  to  the  Rayleigh  distribution.  We  observed 
the  presence  of  a  few  very  intensive  harmonics  (so- 
called  oligarchs,  cf.  ini),  which  did  not  obey  the 
Rayleigh  statistics,  but  their  contribution  to  the 
total  balance  of  the  wave  action  is  small  (no  more 
than  5%).  This  means  that  we  almost  overcame 
negative  effects  caused  by  the  finite  size  of  our  sys¬ 
tem  (see  [H  0,  H3),  and  that  our  grid  is  fine 
enough. 


2.  We  observed  the  formation  of  the  Zakharov- 
Filonenko  spectral  tail  in  the  energy  spectrum  |ryw|2 
(see  Fig.  [3). 


to 


FIG.  3:  The  time-averaged  spectrum  of  a  point  on  the  surface 
nw  =  {\i]u |2)  in  the  double  logarithmic  scale.  The  tail  of  the 
distribution  fits  to  the  Zakharov-Filonenko  spectrum  u>~4. 


At  the  same  time,  we  observed  a  manifestation  of  strong- 
turbulent  effects.  They  are  manifested  by  the  formation 
of  “fat  tails”  on  the  PDF  for  surface  elevations  and  espe¬ 
cially  for  its  gradients  (see  Fig[4j)  The  presence  of  these 


FIG.  4:  PDF  for  (Vrfiy  at  the  moment  of  the  maximum  sur¬ 
face  roughness,  t  ~  147o- 

tails  indicates  that  the  surface  has  a  tendency  to  become 
rough  and  to  produce  white-capping.  In  our  model,  wave¬ 
breaking  is  arrested  by  the  strong  pseudo- viscosity. 


STATISTICAL  EXPERIMENTS 


In  the  second  experiment,  we  solved  the  Hasselmann 
kinetic  equation  for  nj:  =  (|a*;|2)  §1 


dnr  , 

'  =  Snl  [n\  +  Sdiss  T 


fc! 


nr  nr  nr  - 
k 1 


Sni[n]  =2ng2  J  ( 

+nknk2nk3  -  nkVk2  ~  nkVk3] 

x8{uk  +cokl  -  Wfc2  -u>k3)  x 
x S  (  k  +  k\  —  Ji2  —  )  d/cid/c2dfc3. 


(5) 


is  the  phe- 
the  wlrite- 


Here  y*,  is  the  pseudo-viscosity  and  Sdiss 
nomenological  dissipation  term  modelling 
capping  process. 

Eq.  (J5|  was  solved  on  the  grid  71  x  36  in  polar  coor¬ 
dinates  on  the  frequency-angle  plane  by  the  Resio- Tracy 
code  27 1 ,  improved  in  [5,  lj|  ■  We  first  performed  the  ex¬ 
periment  with  Sdiss  =  0.  We  observed  good  qualitative 
coincidence  with  the  dynamical  experiment.  We  observed 
a  downshift  of  the  spectral  peak,  angular  spreading  and 
the  formation  of  uC4  spectral  tails.  But  the  quantitative 
agreement  of  the  experiments  was  not  good:  it  was  clear 
that  the  inclusion  of  some  phenomenological  dissipation 
is  necessary. 

We  examined  the  standard  from  of  Sdiss  used  in  in 
the  industrial  operational  models  of  wave  forecasting  — 
WAM  Cycle  3  and  WAM  Cycle  4  (hereafter  WAM3  and 
WAM4)'0: 


Sdiss  —  Cds&Y  (  (1  ~  8)  +  8 


(6) 


where  k  and  u>  are  the  wave  number  and  the  frequency, 
tilde  denotes  the  mean  value;  Cds,  8  and  p  are  tun¬ 
able  coefficients;  S  =  k\fH  is  the  overall  steepness; 
Spm  =  (3.02  x  10-3)1/2  is  the  value  of  S  for  the  Pierson- 
Moscowitz  spectrum  (note  that  the  characteristic  steep¬ 
ness  is  p  ~  \/2 S).  It  is  worth  noting  that  according  to 
0  ,  the  theoretical  value  of  the  steepness  for  the  Pierson- 
Moscovitz  spectrum  is  Spm  —  (4.57  x  10~3)12,  which 
gives  us  p  ~  0.095. 

The  values  of  tunable  coefficients  in  the  WAM3  case 
are: 


Cds  =  2.36  x  10“5,  8  =  0,  p  =  4  (7) 

and  in  the  W AM 4  case  are: 


Cds  =  4.10  x  10"5,  8  =  0.5,  p  =  4  (8) 

The  evolution  of  the  total  wave  action  is  presented  on 
Fig.  [5]  One  can  see  that  in  the  long  run,  the  models 
WAM3  and  WAM4  overestimate  white-capping  dissipa¬ 
tion.  To  achieve  better  agreement  of  both  experiments, 
we  used  the  following  form  of  the  dissipative  term: 

Cds  =  1-00  x  10~6,  <5  =  0,  p  =  12  (9) 
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The  total  wave  action  curve  corresponding  to  this  new 
dissipation  term  is  shown  on  Fig.[5]by  the  thick  solid  line 
and  displays  excellent  correspondence  with  the  dynami¬ 
cal  model. 


o 


Initial  Wave  Periods 
1000  2000 


3000 


FIG.  5:  Total  wave  action  as  a  function  of  time.  The  solid  line 
corresponds  to  the  dynamical  equations,  the  dashed-dotted 
line  —  to  the  kinetic  equation  with  artificial  viscosity,  the 
dashed  line  —  to  the  kinetic  equation  with  the  WAM3  damp¬ 
ing  term,  the  dotted  line  —  to  the  kinetic  equation  with  the 
WAM4  damping  term,  and  the  thick  solid  line  —  to  the  ki¬ 
netic  equation  with  the  new  damping  term, 


CONCLUSION 


Our  experiments  can  be  interpreted  as  a  confirmation 
of  the  theory  of  weak  turbulence.  However,  even  at  mod¬ 
erate  values  of  the  parameter  of  the  nonlinearity  /z,  the 
strongly  nonlinear  effects  of  white-capping  are  essential. 
They  manifest  themselves  as  fat  tails  of  the  PDF  and 
lead  to  additional  dissipation  of  wave  energy.  This  dis¬ 
sipation  demonstrates  a  very  strong  dependence  on  the 
steepness.  At  steepness  /z  =  0.176  they  dominate,  at 
steepness  fi  =  0.09  they  are  negligibly  small.  The  re¬ 
sults  of  our  experiments  are  in  good  qualitative  agree¬ 
ment  with  the  field  experiment  of  Banner,  Babanin  and 
Young  [l2| ,  who  found  that  wave- breaking  is  a  threshold 
effect,  similar  to  a  second-order  phase  transition. 

We  stress  that  the  dependence  (HI)  is  much  sharper  than 
it  is  usually  stated.  So  far,  the  sharpest  dependence  p  =  5 
was  given  by  Donelan  [28| .  We  can  guess  that  the  real  de¬ 
pendence  of  Sdiss  on  ^  is  even  stronger,  and  that  the  on¬ 
set  of  the  wave  breaking  is  a  threshold-type  phenomenon 
like  a  second-order  phase  transition. 
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